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Abstract 



(— I ! We analyze the latest published solar neutrino data assuming an arbitrary 



neutrino oscillation/conversion mechanism suppresses the electron neutrino 



^ I flux from the Sun independent of energy. For oscillations/transitions into 

active (sterile) neutrinos such mechanisms are ruled out at 99.96 (99.9997) 
% C.L. assuming the standard solar model by Bahcall and Pinnsoneault '95 
correctly predicts all solar neutrino fluxes within their estimated uncertainties. 
Even if one allows for and ^Be solar neutrino fluxes that are vastly different 
from the ones in contemporary standard solar models these mechanisms are 
strongly disfavored by the data. 



*Also at: Institute of Nuclear Research and Nuclear Energy, Bulgarian Academy of Sciences, 1784 
Sofia, Bulgaria. 





I. INTRODUCTION 



The present solar neutrino data are in significant disagreement with the most 

elaborate contemporary standard solar models (SSM) Theoretical uncertainties and 

experimental errors cannot account for this discrepancy and the data indicate that the solar 
Vf. flux is reduced with respect to its predicted value in the SSM. Assuming the electron neu- 
trinos (z/g) produced in the Sun and later detected in the solar neutrino detectors on Earth, 
travel unaltered between the points of production and detection, and using the luminosity 
constraint one finds that the inferred from the data reduction of the intermediate energy ^Be 
neutrinos is much stronger than the reduction of both the low-energy neutrinos, the flux 
of which actually appears to be increased by 10 % at the best fit point, and the high-energy 
®B neutrinos, the flux of which turns out to be lower by a factor of ^ 0.3 

Presently the only known way to explain the discrepancy between the predictions of the 
solar models and the experimental results is to resort to neutrino flavor conversion mech- 
anisms (MSW vacuum oscillations [|10|, etc.) which can provide an energy dependent 
suppression of the known solar z/g fluxes (pp, ^Be, ®B, pep, ^^N and ^^O). An important 
role in the explanation plays also the z/^ (and/or v^) neutral current contribution to the 
neutrino-electron scattering rate, present in the case of solar v^. oscillations/transitions into 
active neutrinos. This contribution increases the ratio of the event rate in the Kamiokande 
detector to the event rate expected from the SSM without oscillations, with respect to the 
analogous ratio for the chlorine experiment, in accordance with the data. 

In a previous publication |]ll| we have shown that neutrino conversion mechanisms which 
predict an energy independent suppression Q of the entire solar neutrino spectrum are ruled 



*Here and in what follows we always have in mind solar neutrino conversion mechanisms which 
lead to an energy independent and equal suppression of all components of the solar v^. flux. For 
brevity we call such mechanisms "energy independent mechanisms". Let us note that, for instance, 
the resonant solar neutrino transitions due to the interplay of neutrino flavour changing and flavour 
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out by the data, assuming the standard solar model [|T^] correctly predicts all neutrino 



fluxes. We also showed that even if the neutrino flux is reduced by one half and the ^Be 
neutrino flux is reduced by 30 % with respect to the corresponding theoretically predicted 
fluxes in the SSM, these mechanisms still provide a poor flt to the data. Examples of such 
mechanisms with energy independent suppression are neutrino spin-precession |]T^ , vacuum 
oscillations with an oscillation length much shorter then the Sun - Earth distance, adiabatic 
MSW transitions for all neutrino flux components, possible only for large mixing angles and 
10~^ eV^ < Am^ < 10"^ eV^ etc. 

Our previous analysis has been published as a part of a larger investigation which was 
centered on somewhat different problems, used older data and by now an outdated reference 
solar model. Here we update our previous results by using the latest solar neutrino data 
summarized in Table. |I[ We extend our analysis by using neutrino fluxes from three recent 
solar models 0-0] which include helium and heavy elements diffusion, an effect that has been 
shown lately to considerably improve the agreement with helioseismological very precise 
measurements of low degree p-modes []T3|,|IB|. We also use recently updated cross-sections 



for the chlorine and gallium detectors ||T^ which improve the accuracy of the predicted event 



rates in these detectors. We explain in detail the assumptions and numerical procedures used 
in the analysis. Our latest results conflrm and strengthen the results and the conclusions 
reached in our previous publication ||TT|. 



II. NUMERICAL PROCEDURE 

Given the published experimental results, a standard solar model which predicts the 
neutrino fluxes, and a mechanism depleting the fluxes of solar by a factor P which 



diagonal (but flavour dependent) neutral current interactions |12| are in the case of massless neutri- 
nos energy independent. Nevertheless, the corresponding transition probabilities can be different 
for the '^B, '^Be, pp, etc. neutrinos (see ref. |12| for details). 
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is independent of energy, one can compute including the theoretical uncertainties and 
experimental errors as described, e.g., in |T^. The theoretical uncertainties are due to 
uncertainties in the input parameters of the standard solar model as well as to uncertainties 
in the detection cross-sections. The correlations between the solar neutrino fluxes which 
arise because one and the same input parameter can affect several neutrino fluxes also have 
to be taken into account. In this calculation we have only one independent parameter, P, 
which is supplied by the relevant mechanism proposed as a solution of the solar neutrino 
problem. We describe the results of this calculation in the next section. 

We next relax the uncertainty intervals, estimated within the SSM, of two of the most 
uncertain solar neutrino fluxes, and ^Be. For three representative values of the beryllium 
neutrino flux, ^^e, equal to 0.7, 1.0 and 1.3 times the flux in the SSM we vary both P 



and the total boron neutrino flux, treating them as free parameters |TT[. This proce- 
dure was adopted in order to study how much model dependent are the strong constraints, 
obtained with the first procedure described above, on the energy independent neutrino 
oscillation/conversion mechanisms. One of the largest uncertainties of the solar model pre- 
dictions is related to the fact that the cross-section for the ''Be(p, 7)^B reaction has never 
been measured at the typical energies at which this reaction is taking place in the Sun. This 
is the primary source of the large existing uncertainty in the predicted boron neutrino flux. 
In the model [§] the estimated one sigma relative uncertainty in the ^B flux is (!tai7) and 
the corresponding uncertainty in the ^Be neutrino flux is (^o o?)- The relative uncertainties 
in the SSM prediction of the CNO neutrino fluxes, ^'^N and ^^O, are even larger than those 
of the boron flux. However, because of their lower energies and relatively low fluxes, the 
CNO neutrinos contribute much less to the event rates in the chlorine and gallium neutrino 
experiments than the ®B and/or ^Be neutrinos and thus are less important in the analyses 
of the solar neutrino data. Therefore in our calculations we keep the CNO neutrino fluxes 
equal to their SSM values. 

The predicted event rates in the operating detectors can be expressed as functions of the 
corresponding event rates as calculated from the SSM, the constant suppression factor, P, 



and the assumed ratios, /„, of each solar neutrino flux, to the same flux in the SSM 0, 
^a,ssM, fa = ^a/^a,ssM, « = PP, ''Be, pep, ^B, ^^N and ^^O. The three equations which 
give these event rates are: 



Qca = P [JbQL + fseQ^: + UpQZ + fpepQ'c^ + QcaiCNO)) , (la) 
QCI = P {IbQEi + fBeQEl + fpepQ'cT + Qci{CNO)) , (lb) 
Ke = fBiOM5PB + 0.155), (Ic) 



or 



Rue = fsPs. (Id) 

The quantities in the lefthand side of the first two equations are the rates of ^^Ge and ^'^Ar 
production in the gallium and chlorine detectors and R^e in Eq- (0) and Eq. ([THD is the 



ratio (measured flux)/(SSM ®B flux) for the Kamiokande (z/ec) scattering experiment. 
Depending on the particular mechanism involved, P can be the solar Uf. survival probability 
for neutrino oscillations, spin-flavor conversion, etc. For neutrino oscillations/transitions into 
active neutrinos (z/^ and/or z/^-), or antineutrinos (z/^ and/or z/^-), Eq. (|T^) is the relevant 
one, whereas Eq. (^) has to be used for oscillations/transitions into sterile neutrinos (z/^). 
Qf is the contribution of source a to the total event rate in the i-th radiochemical detector 
{i = Ga, CI) as predicted in the SSM. The coefficients /„ are written explicitly only for those 
fluxes for which we have assumed that they might be different from the ones predicted in 
the SSM. 

For fse = 0.7 (1.3) we chose to satisfy the luminosity constraint including small correc- 
tions for thermal motion of the particles in the plasma and inequalities between the neutrino 
fluxes, which follow from the sequencial process of nuclear fusion in the Sun |jT9| [], by in- 



creasing (decreasing) the pp neutrino flux by ~ 2 %. We always kept the ratio of the pep 
to pp flux equal to the one in the SSM. The prediction for this ratio in the different solar 



^For earlier discussions of the luminosity constraint see ref. |2C| 
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models is remarkably stable: it doesn't vary by more than one percent. The numerical coef- 
ficients in Eq. (|Ic|) have been computed using the neutrino-electron scattering cross-sections 
with radiative corrections The ratio (0.155) of the (integrated over the ®B spectrum) 



scattering cross-sections cr(z/^(z/T-)e)/o"(z/ee) has been computed using the published trigger 
efficiency functions and energy resolution of the detector. We have taken into account also 
the known evolution of the trigger in Kamiokande . During the first 449 days the detectors 
has been operated with a threshold of 9.3 MeV, later on the threshold has been reduced 
to 7.5 MeV (794 days) and 7.0 MeV (836 days). A slight change in the trigger efficiency 
function during the Kamiokande-II stage due to increased light collection efficiency has been 
neglected. This effect might change the time-averaged ratio at most by a few percent, which 
will not change significantly any of our results. 
Next we computed the defined as: 

The index i denotes the four operating experiments. The statistical and systematic errors 
from each experiment have been combined quadratically to find the total experimental error, 
CTexp{i)- The theoretical error is due to the uncertainties in the solar detection cross- 



sections for the two radiochemical experiments [|T^ and has been calculated as described in 
18| . The theoretical uncertainty for the Kamiokande experiment has been assumed to be 
zero. Here we neglected the theoretical uncertainties related to the input parameters in the 
solar model but instead varied two of the most uncertain fiuxes, $b and $_Be, within much 
broader intervals than their estimated 3ct uncertainty ranges in the SSM. This approach 
allows us to do the same analysis with three independently developed standard solar models 

III. RESULTS 



Assuming the standard solar model ||^ correctly predicts, within the estimated theoretical 
uncertainties, the fiux from each solar neutrino source, the energy independent mechanisms 
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provide a very poor fit to the data. For oscillations/transitions into active neutrinos the 
minimum is 13.0 (1 d.f.) which rules out this scenario at 99.96 % CL. The constraint 
on oscillations/transitions into sterile neutrinos are even stronger. The minimum is 23.5 
which rules out this hypothesis at 99.9997 % CL. The reason for the difference between 
the CL. at which oscillations/transitions into active and into sterile neutrinos are ruled 
out is the absence in the latter case of neutral current contribution to the neutrino-electron 
scattering rate which exacerbates the difficulty in describing simultaneously the Homestake 
and Kamiokande data. 

The constraints given above depend on the solar model used in the analysis. Of the 
three recent solar models only Ref. gives a complete description of the estimated 
uncertainties in the neutrino fluxes related to the inaccurately known input parameters. In 
order to study this dependence on the solar model, as explained in the previous section, we 
disregard the solar model uncertainties of the most uncertain of the major solar neutrino 
fluxes, $B and ^Be- For j^e we chose three representative values, jse = 0.7, 1.0 and 1.3. In 
each of these three separate cases js as well as P are treated as free parameters. For each of 
the three solar models |^-|^ we computed the (Eq- 0) for a large number of pairs (/s, P) 
from the following intervals: < < 2 and < P < 1, and found the minimum, Xmm) 
for each of the three representative values of jse- The value of Xmm determine the CL. at 
which the tested hypotheses are ruled out by the data. The results are summarized in Table 
|I|for the solar models and for oscillations/transitions into active as well as into sterile 
neutrinos. 

For /Be = 1 the energy independent solar v^, suppression mechanisms are ruled out at 
(96.5)% CL. (x^ = 6.8 (2 d.f.)) for all three solar models studied here. Even if we relax the 
constraint on the ''Be flux, the energy independent mechanisms still give a poor fit to the 
data. The smallest minimum Xmm — 5-9 for oscillations/transitions into active neutrinos is 
achieved for j^e = 0.7. This rules out any energy independent mechanism as a solution of 
the solar neutrino problem at 94.5 % CL. (2 d.f.). Higher values of jse are ruled out at even 
higher (~ 97.5 % for j^e = 1-3) CL. The constraints on neutrino oscillations/transitions 
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into sterile neutrinos are somewhat stronger, e.g., for /bs = 0.7 the Xmm is 8.2. This rules 
out such mechanisms at (at least) 98 % CL. 

The week dependence of the results on the assumed beryllium neutrino flux shows that 
the residual theoretical uncertainties, which we have neglected here and which are related to 
all other neutrino fluxes except and ''Be, cannot significantly change these results. Those 
uncertainties are either considerably smaller {pp and pep neutrinos) than the uncertainties 
in the ^B and ''Be neutrino fluxes, which have been allowed to vary outside their estimated 
3a limits, or are related to sources (CNO), which give minor contributions to the event 
rates in the operating experiments. In several test cases we have verified this by a complete 
calculation where for the model 0, in addition to varying fs and P, we took into account 
the residual uncertainties in the CNO, pep and pp neutrino fluxes. 

In our analysis we did not include solar models without diffusion, e.g., [^]. These models 



are strongly disfavored by helioseismological data [0,0 although the predicted neutrino 
fluxes are only slightly lower than in the more advanced models included in Table We 
also did not include results for the model in ref. [^. The results for this model, which differs 
from the models |§-0] mainly by the choice of the astrophysical factors, give somewhat lower 
Xmin cases changing the CL. by less than 0.5a. 

The constraints so obtained on the energy independent solar suppression mechanisms 
depend only weakly on the particular solar model used as a reference. In addition, we 
want to emphasize that at present there is no compelling reason to expect such strong 
deviations, as the ones considered here, from the predicted fluxes in standard solar models. 
Our goal was to show that even ad hoc changes in the solar neutrino fluxes cannot avert 
our main result, namely that the neutrino oscillation/conversion mechanisms which predict 
energy independent suppression are strongly disfavored by the data from the operating solar 
neutrino experiments. 

In order to test the significance of the data from each solar neutrino experiment for 
the poor fit provided by energy independent neutrino oscillations/transition mechanisms we 
have computed the for each of the three possible pairs of experiments, namely (ClAr + 



(z/e)), ((i^e) + GaGe) and (GaGe + ClAr). When only P is allowed to vary and the neutrino 
fluxes are assumed to be the ones from the SSM in Ref. 0, the minimum for the three 
cases, assuming oscillations/transitions into active neutrinos, are: 



ClAr + (i/e) : xLn = 0.07; P = 0.28, (3a) 

(z/e) + GaGe : xLn = 2-22; P = 0.50, (3b) 

GaGe + ClAr : xLn = 10.7; P = 0.46. (3c) 

For oscillations/transitions into sterile neutrinos we obtain: 

ClAr + (z/e) :xLn = 1-15; P = 0.34, (4a) 

(z/e) + GaGe : xLn = 0-79; P = 0.50, (4b) 

GaGe + ClAr : xLn = 10.7; P = 0.46. (4c) 



From these results it is evident that (both for oscillations/transitions into active and sterile 
neutrinos) when only (GaGe + ClAr) experiments are considered the energy independent 
solutions are strongly disfavored by the data whereas both (ClAr + (z/e)) and (GaGe + 
(z/e)) data can be fltted by an energy independent suppression. 

When both P and the boron neutrino flux are allowed to vary for any pair of experimental 
results there is always an acceptably small Xmin- The values of P and fs for which is 
minimum are slightly different for the different solar models. However, the assumption that 
both the SSM and at least one of the solar neutrino experiments are wrong, seems too 
extreme and hardly justiflable. 

Some energy independent mechanisms predict a flxed value of the survival probabil- 
ity which cannot be varied. Examples of such mechanisms are neutrino oscillations with 
maximal mixing between two (P = 1/2) or three (P = 1/3) neutrinos |2^, as well as the 
mechanism discussed in ref. (P = 5/9). The results in Table |ITT| show that all three 
values of the survival probability, P = 1/3, 1/2, 5/9, are ruled out by the data both for 
oscillations/transitions into active and into sterile neutrinos. The Xmin values found corre- 
spond respectively to 99.4, 96.5, 94.5 % C.L. and to 99.96, 99.4, 98.4 % C.L. for the two 
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types of oscillations/transitions. Not only are the minimum too high for all three models 
considered, but the values of the boron and especially of the beryllium neutrino fluxes, for 
which the CL. is lowest, are far outside the acceptable uncertainty ranges for these fluxes. 

IV. CONCLUSIONS 

Mechanisms proposed as solutions of the solar neutrino problem and predicting an en- 
ergy independent suppression of the solar neutrino spectrum as a result of neutrino oscil- 
lations/transitions into active (sterile) neutrinos or antineutrinos are ruled out by the data 
from the four operating experiments at 99.96 (99.9997) % C.L., if the standard solar 
model 1^ correctly predicts (within the estimated uncertainties) the solar neutrino fluxes 
from the different sources. Even relaxing the constraints on the ''Be neutrino flux by assum- 
ing that it differs by 30 % from it's values in three recent solar models with heavy element 
diffusion, and treating the neutrino flux as a free parameter, still provides a poor fit to 
the data for any value of the energy independent suppression factor. Suppression factors of 
1/3, 1/2, and 5/9 0,|5[ are ruled out by the data at 99.4, 96.5, 94.5 % CL. (99.96, 99.4, 
98.4 % CL.) for oscillations/transitions into active (sterile) neutrinos even if one treats both 
the ^B and ^Be fluxes as free parameters in the analysis. Not only are the minimum too 
high, but the lowest Xmm values are reached for boron and especially beryllium neutrino 
fluxes which lie far outside the acceptable uncertainty ranges for these fluxes. 
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TABLES 

TABLE L Solar neutrino data used in the analysis and the theoretical predictions from the 
SSM m. The units are SNU for the event rates in the radiochemical detectors and cm~^s~^ for the 
flux measured by Kamiokande. For both the chlorine and gallium experiments the theoretically 



predicted event rates have been updated following the results in Ref. |17|. 



Experiment 



Data 



Theory 



Units 



Chlorine 

Kamiokande 

GALLEX 

SAGE 



2.56 ±0.16 ±0.14 
2.80 ±0.19 ±0.33 



9.5+1-2 



6.62 



1.4 



hO.93 
-1.12 



SNU 



10^ cm-^s-i 



69.7 ± 6.7" 



-3.9 
-4.5 



136.8 



-7 



72.0 



-12-1-5 
-10-7 



136.8 



+8 
-7 



SNU 
SNU 
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TABLE II. Minimum ^ for the energy independent mechanisms including oscilla- 
tions/transitions into active or sterile neutrinos. The results are given for three reference solar 
models fSj-Q. The ''Be neutrino flux has been chosen: a) equal to the one in the corresponding 
solar model, b) 30 % higher and c) 30 % lower. The values of /s, the ratio of the boron neutrino 
flux to the predicted ^B flux in the relevant SSM and P, the energy independent suppression factor, 
for which Xmm achieved in each case, are given also. 









active 






sterile 




solar model 






fB 


P 




fB 


P 




0.7 


6.6 


0.47 


0.55 


9.2 


0.47 


0.56 


BP'95 


1.0 


7.5 


0.51 


0.49 


10.7 


0.49 


0.51 




1.3 


8.2 


0.54 


0.45 


12.0 


0.51 


0.47 




0.7 


5.9 


0.61 


0.57 


8.2 


0.60 


0.58 


CDFLR 


1.0 


6.8 


0.64 


0.52 


9.6 


0.62 


0.54 




1.3 


7.5 


0.67 


0.48 


10.8 


0.64 


0.50 




0.7 


6.7 


0.50 


0.56 


9.3 


0.49 


0.57 


RVCD 


1.0 


7.6 


0.54 


0.50 


10.9 


0.52 


0.52 




1.3 


8.4 


0.57 


0.46 


12.3 


0.54 


0.48 
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TABLE III. Minimum and the corresponding values of and f^e for three characteristic 
values of P (the solar Vf, survival probability) in specific models predicting energy independent 
suppression of the solar neutrino spectrum. The results are given for three solar models (Bj-Q. In 
each case the upper row is for oscillations/transitions into active neutrinos, whereas the lower row 
is for oscillations transitions into sterile neutrinos. 







P = 1/3 






P = 1/2 






P = 5/9 




solar model 




/b 






fB 






fB 


fBe 


BP'95 


10.3 


0.71 


2.0 


7.0 


0.53 


0.75 


6.0 


0.50 


0.44 




16.2 


0.74 


2.0 


10.1 


0.55 


0.67 


8.3 


0.52 


0.36 


CDFLR 


11.1 


0.97 


2.0 


6.8 


0.69 


0.91 


5.8 


0.65 


0.55 




16.5 


1.01 


2.0 


9.8 


0.72 


0.81 


8.1 


0.67 


0.48 


RVCD 


11.1 


0.78 


2.0 


7.4 


0.57 


0.78 


6.3 


0.53 


0.48 




17.2 


0.81 


2.0 


10.6 


0.59 


0.70 


8.8 


0.55 


0.40 
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